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In spacetimes of dimension greater than four it is natural to consider higher order (in R) corrections 
to the Einstein equations. In this letter generalized Israel junction conditions for a membrane in 
such a theory are derived. This is achieved by generalising the Gibbons-Hawking boundary term. 
The junction conditions are applied to simple brane world models, and are compared to the many 
contradictory results in the literature. 
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It is well known that for a boundary-less spacetime, 
the Einstein equations can be derived from the Einstein- 
Hilbert action. This is no longer true if the spacetime 
has a boundary. The problem is resolved by adding a 
Gibbons-Hawking boundary term to the action . Vary- 
ing the action then gives the correct held equations, as 
well as boundary conditions at the edge of the spacetime. 

A slight variation of this is to consider an infinitely 
thin (D — 2)-brane in a D dimensional spacetime. Since 
spacetime is split into two, the brane can then be treated 
as the boundary of each half of the spacetime. Varying 
the Gibbons-Hawking term then gives the Israel junction 
conditions on the membrane |2j. These junction condi- 
tions have recently received a great deal of attention due 
to their use in the study of 'brane- worlds' (see for exam- 
ple (|, ||, H, H 0]). In a brane- world scenario our universe 
is modeled by a 3-brane embedded in a five-dimensional 
'bulk' spacetime. In the simplest cases, all the standard 
model fields are confined to the brane, while only gravity 
propagates in the bulk. The brane is usually taken to be 
of zero thickness, and so the Israel junction conditions 
can be used to relate the bulk dynamics to what we, on 
the brane, observe. 

The bulk gravitational field equations are usually as- 
sumed to be the five-dimensional Einstein equations. 
However in spaces of dimension greater than four it is 
natural to consider additional higher order curvature 
terms § g |fij 0. 

In general relativity, the vacuum field equations are 
taken to be a linear combination of the Einstein tensor 
and the metric. This choice is motivated by the fact 
that it is the most general combination of tensors which 
is (a) symmetric, (b) depends only on the metric and 
its first two derivatives, (c) is divergence free and (d) is 
linear in second derivatives of the metric. In fact, in four 
dimensions, the fourth condition is superfluous since it 
is implied by the other three [g. In five dimensions the 
second order Lovelock tensor 

H a b = RRab — 2R ac R C b — 2R cd R acb d + Ra Cde Rbcde 

— -^9ab {R 2 — 4:R c dR cd + R cdes Redes) (1) 

also satisfies the above conditions. Thus the most general 
choice of gravitational vacuum field equations in five di- 
mensions is a linear combination of g abl G ab and H ab . In 



the absence of any experimental evidence to the contrary, 
all three terms should be included. 

A further motivation for higher order curvature terms 
is that they also appear in the low energy effective field 
equations arising from most string theories. Since brane 
worlds are motivated by string theories, it is particularly 
natural to include the extra terms in the five-dimensional 
field equations. 

Like the standard Einstein equations, these higher or- 
der generalisations can also be derived from an action. 
The tensor H ab can be obtained from an action contain- 
ing the Gauss-Bonnet term 

£gb = R 2 - 4R ab R ab + R abcd R ab ed ■ (2) 
Consider the action 

S M = i^2 I d °xV^9 {R-2A + <xC G b} 

for a D dimensional manifold M.. In a string theory 
context we would have k~ 2 — M®~ 2 and a oc M~ 2 , 
where M* is the string mass scale. 

Let us suppose, as would be the case in a co-dimension 
one brane world scenario, that M. is split into two parts 
by a hypersurface S, whose two sides will be denoted £±. 
Their normals, n a , will be taken to point away from the 
surface and into the adjacent space. 

Varying the action (g) with respect to the metric gives 



(3) 



SSm 



[ d D x^g~8g ab (G ab + Ag ab + 2aH ab ) 

flY ]b + 2aP abcd )v d 6g bc (4) 




where h ab = g ab — n a n b is the induced metric on S, and 

Pabcd = Rabcd + ^Rb[c9d\a ~ 2i? a [cffd]b + Rg a [c9d]b (5) 

is the divergence free part of the Riemann tensor. 

Expression (^) contains normal derivatives of the met- 
ric variation. As with the Einstein-Hilbert action, we 
must add a boundary term in order to cancel them. For 
an action with a Gauss-Bonnet term ^ , the appropriate 
term is fttOl 



2± 



d D - x x^h [K+ 2a{J - 2G ab K ab } 



(6) 
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where K is the trace of the extrinsic curvature, defined 
by K ab = h c a V c n b , and J is the trace of 



1 



J ab =-[2KK ac K % + K cd K^K ab 

- 2K ac K cd K db - K 2 K ab ) , (7) 

Throughout this paper I will denote tensors associated 
with the induced metric by a circumflex accent, so G ab 
is the (D — 1) dimensional Einstein tensor on E corre- 
sponding to h ab . 

Varying the action Sm + now gives an expression 
which does not contain normal derivatives of 5g ab . If we 
also include a matter contribution to the action 

Smat = - / d?Xy/=g£M- f dP- l xy/=hCW , (8) 

J M JT, 

then the variation of the total action S — Sm + S% + >S ma t 
gives 



G ab + 2aH ab + Ag ab 



(9) 



in A4, and 



2(K ab -Kh ab )+4a(3J ab - Jh ab + 2P acdb K cd ) = -n 2 S ab 

(10) 

on S, with (X) = [X(S+) + X(£_)]/2 denoting the 
average of a quantity over the two sides (£±) of the 
hypersurface. The two energy-momentum tensors are 
defined by T ab = 26£ { £ i) /Sg ab - g ab L { ^ ] and S ab = 
25£ < £ ] /6h ab -h ab £% ) . 

With the aid of the Gauss-Codazzi equations (|2~4|- 
p6|) below, we obtain the energy-momentum conservation 
equation on the hypersurface 

D a S ab = -2 (n a (G ac + 2aH ac )h%) = -2k 2 (n a T ac h c b ) , 

(11) 

with D a denoting the covariant derivative corresponding 
to h ab . This is very similar to the corresponding result 
in the standard brane world models 

Many previous papers have tried to derive the junction 
conditions by treating the hypersurface £ as a 5-function 
contribution to the to £^ V< \ as in the original brane cos- 
mology papers f|. In this case there is some ambiguity 
as to the correct definition of H ab on £, which has led 
to the suggestion that the junction conditions must de- 
pend on the thickness of the brane |L2| . This would be 
true for a general combination of second order curvature 
terms, whose action would contain third (or higher) or- 
der derivatives of the metric. However, this is not true 
for the Gauss-Bonnet combination (Q) since it has been 
specifically chosen not to contain such derivatives. 

By allowing n c d c g ab to be discontinuous at the brane, 
and treating (n c d c ) 2 g ab as a 5-function, junction condi- 
tions which are independent of the brane thickness can 
be found jD| [H], [l5j [l(| |l^|. However, care must be 
taken to regularise the <5-function correctly |l3|, [l8| . This 



was done in refs. [^3[ Q and the resulting junction 
conditions agree with those in is letter. Refs. @, do 
not regularise the J-function appropriately, and obtain 
incorrect results. 

We will now use the above results to derive the Fried- 
mann equation for a cosmological brane world model. 
As has been shown by Kraus and Ida H (for the stan- 
dard five dimensional Einstein equations), it is possible 
to obtain a cosmological generalisation of the Randall- 
Sundrum || model by considering a static bulk space- 
time, and allowing the brane to have a time dependent 
position in the bulk. 

The bulk metric can be written in the form 



ds 2 



-h(r)dT 2 



dr 2 
h(r) 



r 2 ^lijdx % dx^ 



(12) 



where fly is the three dimensional metric of space with 



constant curvature k = —1, 0, 1. For C 
equations are solved by [O 



(M) 



= 0, the field 



r 

4a 



4 u 
1 + -aA + 8a^j 
3 r 4 



(13) 



with /i being an arbitrary constant. In the a — * limit, 
fj, is equal to the black hole mass. 

We define the position of the brane as r — a(r) and 
T = t(r), which is parameterised by the proper time on 
the brane t. The induced metric is 



ds 2 



-dr 2 + a(r) 2 flijdx i dx^ , 



(14) 



the tangent vector of the brane is u a — (T, 0, 0, 0, f) and 
n a = (— f, 0, 0, 0, T). Normalisation of n a implies 



h 2 T 2 



= -h 



(15) 



We take the brane matter to be a perfect fluid, so S ab — 
(p + P) u a u b + ph ab - The (uu) component of ( |l0| ) is then 



1 + -aH 2 + 4a- 



(hf) 4 (h 2 T) p 



-a 



(16) 



For simplicity let us assume Z2 symmetry across the 
brane. Squaring ( |l6| ) and simplifying with ( |l5| ) gives a 
cubic equation for H 2 . This has the real solution 



where 



2 k c+ + c_-2 



1 + | aA + 8a ?) 8/a + £ T ± Vf 



(17) 



2/3 



(18) 

This Friedmann equation agrees with rcf. |l3[ , but not 
ref. jl(| [l7| . The only difference between these two con- 
flicting results is one factor of 3, but surprisingly this 
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gives a substantially different Friedmann equation. For 
p ^constant, eq. jl7| ) also agrees with ref. Ref. 
uses different boundary terms [^0) to obtain a different 
Friedmann equation. However the terms used do not give 
a consistent action, and so the result is incorrect. 

As in the usual brane cosmology ||, the standard 
Friedmann equation can be recovered at late time (large 
a) by splitting p into a cosmological constant and a mat- 
ter part Jl3[ |. 

The boundary terms (^) are easily generalised to ac- 
tions where other fields couple to the curvature tensors. 
Variation of the action 

S=A/ d D x^^(x p )R-2A + a^{x p )C GB ) 
2^ J M 

-\ [ d D - 1 xy/^h(^(x^)K 

+ 2a^f(x fi ){ J - 2G ab K ab }) + S mat (19) 

produces the field equations 

<S>G ab - V Q V b $ + 5ah V 2 $ + Ag ab 

+ 2a^H ab + 4aP eacb V e V c * = n 2 T ab (20) 

and the boundary conditions 

{<$>{K ab - Kh ab ) - h ab n e d e $) 

+ 2a(3* J ab - VJh ab + 2^P acdb K cd ) 
+ 2a({2G ab + 2K ea K\ - 2KK ab + 

h ab [K 2 - K cd K cd ]}n e d e V) 

+ K c[a h d]b - Kh a[c h b]d )D c D d V) 

= -\^ 2 S ab . (21) 

All the problematic boundary terms, like those appear- 
ing in ([|), cancel out. In a string theory context and 
\E' would typically be functions of the dilaton or moduli 
fields. 

If we consider a conformally flat spacetime of the form 
ds 2 = -e 2A{y \dT 2 + dxUxi) + dy 2 (22) 

and take $, ^ and Cm = A/k 2 to be functions of y only, 
the generalised Israel junction conditions reduce to 

2 (3<M' + - 4aA' 2 (A'^ + 3*')) = -A (23) 

This agrees with other results in the literature ]15|. 



Appendix 

To prove equation (|Io|), we first use the Gauss-Codazzi 
equations 

R P qrsh p a h q b h r c h s d = R abc d + K bc K ad - K ac K bd (24) 
n a R aqrs h\h r c h s d = D d K bc - D c K bd (25) 

Rm = R a qcp h c a h\h p d + KK bd - K bc K c d (26) 

and contractions of them to expand n a P abcd in terms 
of K ab , n a and quantities associated with the induced 
metric: 

n a P abcd = 2D [d K c]b + 2D e K e [d h c]b + 2h b[d D c] K 

+ 2G b[c n d] + 2(K b e - Kh b e )K e[c n d] 
+ (K 2 - K ae K ae )h b[c n d] . (27) 

The covariant derivative D is defined by D a X bc ... = 
h ah p b h q c • • • V m Apq... . 

To find the variation of Ss (^|) with respect to g ab , 
we first note that the normalisation of n a implies 5n a = 
^n a n c n d 5g cd . Thus, after a little algebra, and using the 
definitions of K ab and D a , we obtain 

h c a h d b 5K cd = n e W [e Sg p]q h p {a h q b) - ^Sg cd K c{a h b)d 

- \D {a (h c b) Sg ce n e ) (28) 

and, from (p6|), 

h c a h d b SR cd = D ia D e (h c b) h d e Sg cd ) 

- ^D e D e (h c a h d b 6g cd ) - ^D a D b (h cd Sg cd ) . (29) 

Now, with the aid of integration by parts and the relation 

Y abe X7 e S 9ab = D e (Y abe 6g ab ) + 2n c 6g cd Y^ e K pe 

+ Sg ab D e Y abe , (30) 

which holds for any tensor Y abe which is orthogonal to 
the normal n c , the variation of the total action S can be 
reduced to eq. (|To|). The more general case (^l]) can be 
dealt with in the same way. 

The energy conservation equation ( [Ti"|) can be proved 
by using D a P acdb = and the Gauss-Codazzi equa- 
tions (p4|-p6|) to express the divergence of (|l0| ) in terms 
of bulk tensors. 



Acknowledgements 



I am grateful to Jihad Mourad and Christos Char- 
mousis for useful and interesting discussions. This work 
was supported by EC network HPRN-CT-2000-00152. 



* Electronic address: Stephen.Davis@th.u-psud.fr 
[1] G. W. Gibbons and S. W. Hawking, Phys. Rev. D 15 
(1977) 2752. 



4 



[2] 

[3] 



W. Israel, Nuovo Cim. B 44 (1966) 1. 
K. Akama, Lect. Notes Phys. 176 



friep-th/OOOlli: 



(1982) 267 



V. A. Rubakov and M. E. Shaposhnikov, Phys. Lett. B 
125 (1983) 136; 

M. Visser, P hys. Lett. B 159 (1985) 22 



|tiep-th/991009: 



[4] 
[5] 
[6] 

[7] 

[8] 
[9] 



[10] 
[11] 



[12] 



E. J. Squires, Phys. Lett. B 167 (1986) 286. 

H. A. Cha mblin and H. S Reall, Nucl. Phys. B 562 

(1999) 133 ]hep-th/9903225 |. 

L. Ra ndall and R. Sun drum, Phys. Rev. Lett. 83 (1999) 
4690 [hep-th/9906064 |. 

P. Binetruy, C. Deffayet and D. L anglois, Nucl. Phys. B 
565 (2000) 269 |hep-th/9905012[ ]; 

P. Binetruy, C. Deffayet, U. E llwanger and D. Langlois, 

Phys. Lett. B 477 (2000) 285 ]hep-th/99102is|]. 

P. Kraus, JHEP 9912 (1999) 01 1 |hep-th/991014s| l; 

D. Ida, JHEP 0009 (2000) 014 |gr-qc/9912002[ . 

D. Lovelock, J. Math. Phys. 12 (1971) 498. 

J. Madore, Phys. Lett. A 110 (1985) 289; Phys. Lett. A 

111 (1985) 283; 

N. Deruelle and J. Madore, Phys. Lett. A 114 (1986) 

185; Phys. Lett. B 186 (1987) 25; 

J. T. Wheeler, Nucl. Phys. B 268 (1986) 737. 

R. C. Myers, Phys. Rev. D 36 (1987) 392. 

D. G. Boulware and S. Deser, Phys. Rev. Lett. 55 (1985) 

2656; 

R. G. Cai, Phys. Rev. D 65 (2002) 084014 



friep-tri/0109133; ]. 

N. Deruelle and T. Dolezel, Phys. Rev. D 62 (2000) 



103502 fer-qc/0004021 1. 

C. Charmousis and J. F. Dufaux, 



[16] 
[17] 

[18] 
[19] 
[20] 



Quant. Grav. friep-th/0202107 |. 

K. A. Meiss ner and M. Olecho wski, Phys. Rev. Lett. 86 



to appear in Class. 



[13 
[14 

[15] I Low and A. Zee, Nucl. Phys. B 585 (2000) 395 



(2001) 3708 [tiep-th/0009122 



[tiep-th/0004124] 



N. E. M avromatos and J Rizos, Phys. Rev. D 62 (2000) 
124004 frep-th/0008074[ ; 

I. P. Neupane, JHEP 0009 (2000 ] 040 |hep-th/0008190 



N. E. Mavromatos and J. Rizos, tiep-th/0205299; 
P. Binetruy, C. Charmousis, S. C. Davis and J. F. Du- 



faux, Phys. Lett. B 544 (2002) 183 fcep-th/0206089 1 
C. Germani a nd C. F. Sopuert a. Phys. Rev 



Eett: 88 



M. Lee, Nucl. Phys. 
ibid. 591 (2000) 587 



(2002) 231101 |hep-th/0202060 | 
J. E. Kim, B. Kyae and H. 
B 582 (2000) 2 96, Erratum: 
Jhep-th/0004005| ]; 
B. Abdesselam and N. Moham medi, Phys. Rev. D 65 
(2002) 084018 |hep-th/0110143| . 

Y. M. Cho, I. P. N eupane and P. S . Wesson, Nucl. Phys. 
B 621 (2002) 388 ]hep-th/0104227t . 

J. E. Lidse y, S. No.jiri and S. D. Odintsov, JHEP 0206 
(2002) 026 ]hep-th/0202~L96 ]. 

A. D. Barvmsk y and S. N. So lodukhin, Nucl. Phys. B 
479 (1996) 305 |gr-qc/9512047} ]; 

M. Cvetic, S. N o.jiri and S. D O dintsov, Nucl. Phys. B 
628 (2002) 295 ]hep-th/0112"045| . 



